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ABSTRACT 

The formation of the first stars is an exciting frontier area in astronomy. Early 
redshifts {z ~ 20) have become observationa lly promising as a result of a recently rec- 
ognized effect ()Tseliakhovich fc Hiratall20f 0[ ) of a supersonic relative velocity between 
the dark matter and gas. This effect produces prominent structure on 100 comoving 
Mpc scales, which ma kes it much rnore fe asible to detect 21-cm fluctuations from the 
epoch of first heating ()Visbal et al.ll2012[ ). We use semi-numerical hybrid methods to 
follow for the first time the joint evolution of the X-ray and Lyman- Werner radiative 
backgrounds, including the effect of the supersonic streaming velocity on the cosmic 
distribution of stars. We incorporate self-consistently the negative feedback on star 
formation induced by the Lyman- Werner radiation, which dissociates molecular hy- 
drogen and thus suppresses gas cooling. We find that the feedback delays the X-ray 
heating transition by a Az ~ 2, but leaves a promisingly large fluctuation signal over 
a broad rcdshift range. The large-scale power spectrum is predicted to reach a max- 
imal signal-to-noise ratio of S/N'^ 3 — 4 at z ~ 18 (for a projected first-generation 
instrument), with S/N> 1 out to 2 ^ 22 — 23. We hope to stimulate additional numer- 
ical simulations as well as observational efforts focused on the epoch prior to cosmic 
reionization. 
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C3 ; 1 INTRODUCTION 

Observations of the redshifted 21-cm line of neutral hydro- 
gen, planned for the next decade, are expected to usher in a 
new era of direct probing of the epoch of first stars. Though 
currently the main observational focus is on the reionization 
epoch, there are instruments hoping to observe the 21-cm 
signal from z ~ 10 - 30, e.g., LEDA (iBurns et al.ll20lj ). 
DAReQ, and the SKA ICarilH et al.ll20oi v 

The formation of the first stars is a relatively clean the- 
oretical problem, as t hey are formed in a metal- free environ- 
ment via H2 cool ing ( Tcgm ark et al.|[l997l : [Machacek et al.l 
I2OOII : I Abel et all [2002). The radiation produced by these 
first radiant objects changed the cosmic landscape dramati- 
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cally (|Madaulll997l '). Three wavelength regimes of this radi- 
ation are most important to consider: the Lyman-a photons 
couple to the 21-cm line at high redshifts (z ^ 30) throug h 
the Wouthuysen-Field effect (|Wouthuvsenll952l : [Fieldlll959l ): 
X-ray photons, produced by stellar remnants, heat the gas; 
and Lyman- Werner (LW) photons dissociate molecular hy- 
drogen, thus produc ing negative feedback on star formation 
(|Haiman et al.|[l997l ) and decreasing the heating rate. 

The radiation spreads out to ~ 100 Mpc around 
each star, where this finite effective horizon arises from 
redsh ift, time delay an d optica l depth effects (|Ahn et al" 
20091 : iMesinger et al.l I2OIII : iHolzbauer fc Furlanetto 



formation progresses, 
up. Fluctuations in 



20 111 : IVisbal et al.l 120121 ). As star 

the radiative backgrounds bui ld 

the radiative backg r ounds (|Barkana fc Loebl l2005l : 
IPritchard fc Furlanetto! I2OO6I ). c aused by the strongl y 
fluctuating distribution of stars (|Barkana fc Loebl |2004| ). 
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couple to the hyperfine transition of neutral hydrogen 
and imprint fluctuations in the redshifted 21-cm signal. 
Whereas Lyman-a coup l ing s at urates at h i gh re dshifts 

the 



Whereas L/yman-a coup l ing s at urates at n i gn re d! 
l|Holzbauer fc Furlanettol I2OIII : IVisbal et all I2OI2I '). 



heating fluctuations couple to the 21-cm brightness tem- 
perature at lower redshifts and thus are more interesting 
in terms of the observational prospects. In particular, 
IVisbal et al.l l|2012l l predict detectible heating fluctuations 
from the first stars at z = 20 with a distinctive signature 
of Baryon Acoustic Oscillations (BAOs) imprinted by 
the supersonic relative (stre aming) velocity between the 
baryons and the dark mat ter (| Tseliakhovich fc Hiratall2010l : 
iDalal. Pen fc SehakllioTol '). 

In this paper we study the signature of the heat- 
ing fl uctuations in cl uding the affect of relative velocity 
as m IVisbal et al.l l|2012l ). and add for the flrst time 
a detailed three-dimensional calculation of the inhomo- 
geneous negative feedback by LW photons. We use the 
same semi-numerical h ybrid methods as in Visbal et al.l 
('2012') ( based in part on i Tseliakhovich fc Hiratal (|201(]| l and 
[iVIcsingo r et al.l l|201ll 'l'l to build a simulation where the stel- 
lar population and the radiative backgrounds evolve simulta- 
neously in time. The basic idea is to linearly evolve a realistic 
sample of the Universe on large scales while using the results 
of numerical simulations and analytical models to add in the 
stars. For a detailed discussion of ou r computational m eth- 
ods, we refer the interested reader to I Visbal et al.l (|2012h . in 
particular to section SI of this paper's Supplementary Infor- 
mation, and references w ithin. We use the sta ndard set of 
cosmological parameters (jKomatsu et al.ll2010h along with 
an assumed star formation efficiency (fraction of gas in star- 
forming halos that turns to stars) of /* — 10%, an X-ray 
photon efficiency of lO^^M^^ based on observed starbursts 
at low redshift as in iMesinger et ahl Hoilj), and LW param- 
eters as explained in the next section. 



2 INCORPORATING THE NEGATIVE 
FEEDBACK 

The formation of the first stars via cooling of molecu- 
lar hydrogen is a highly non-linear pr ocess that can b e 
mimicked by numeri cal simulations, e.g., lAbel et all (|2002h : 
iBromm et ahl l|2002l ). However, numerical simulations in 
which primordial stars are created usually do not consider 
the potentially fatal effect of the LW background on this 
process. The negative feedback of the LW background on 
star formation has been tested in the l i mited case of a fixed 
intensity Jr.w (iMachacek et all I2001I : IWise fc Abel l2007l : 
lO'Shea fc Normanll2008l )~ The feedback boosts the minimal 
cooling mass, Mcooi, i.e., the mass of the lightest halo in 
which stars can form, with the results of these simulations 
well-described by the relation 

Meooi (J21, 2) = Mcooi,o(z) X [l -f- 6.96 (47r Jai)"-''^] , (1) 

where J21 — JLw/(10^'^^erg s~^cm~^Hz~^sr~^) in terms 
of the LW intensity Jlw; another common notation is the 
LW flux -Flw = 47rJLw. Here Mcooi,o(2:) is the value of the 
minimum cooling mass in the standard case with no LW 
background. 

This result is incomplete for two reasons. One is that 
it does not account for the relative velocity Ubc, which 



has a strong impact on the primordial star formation by 
(amo ng other things) bo osting the minimum cooling mass 
(Gre ff et al.ll201ll : IStacAr et al. 2011: Fi alkov et al. 2011,). To 
account for the velocity, we change Mcooi.o(2:) in Eg. ([T]) to 
A/cool, o(-2, ''be), using the fit we developed in iFialkov et al.l 
(201]]) to the streaming- velocity simulations (here we use 
their fit designed for Adaptive Mesh Refinement (AMR) sim- 
ulations as equation (O was the result of a fit to an AMR 
simulation). Thus, we combine two separate physical phe- 
nomena, i.e., the relative motion and the LW fiux, assum- 
ing that they each have a fixed multiplicative effect on the 
minimum cooling mass. This simple ansatz for the depen- 
dence of Mcooi on the two parameters, Vhc and J21, should 
be checked by detailed numerical simulation, which we hope 
to stimulate with this work. We choose this multiplicative 
ansatz since both effects have an independent effect on the 
minimal cooling mass. Therefore they likely reinforce each 
other when they are both present. We note that while we 
include the Vhc and LW effects separately in Eq. (1), our 
results do account for the strong correlation between the ve- 
locity and the LW flux (due to the effect of the velocity on 
star formation). 

The second incompleteness of Eq. U]) is its validity only 
in the case of a fixed background intensity during the forma- 
tion of the halo, whereas in reality the LW intensity is ex- 
pected to rise exponentially with time at high redshifts (e.g ., 
see the Supplementary Information in I Visbal et al.l l|2012l )). 
Treating the intensity as fixed at its final value would greatly 
overestimate the strength of the feedback, since the cool- 
ing and collapse involved in star formation should respond 
with a delay to a drop in the amount of H2. For instance, if 
the halo core has already cooled and is collapsing to a star, 
changing the LW fiux may not stop or reverse the collapse 
at all, and certainly not immediately. Another indication for 
the gradual process involved is t hat the simulation res ults 
can be approximately matched (|Machacek et al.ll200ll ) by 
comparing the cooling time in halo cores to the Hubble time 
(which is a relatively long timescale). Though the relation 
in Eq. ((T)) is the best currently available, more elaborate nu- 
merical simulations, which we again hope to stimulate, are 
needed in order to find a more realistic dependence. We over- 
come this limitation by using the above relation not with the 
final value of J21 at formation, but with the value at a mean, 
characteristic time within the halo formation process. 

The idea of looking at the flux at times well before 
virialization is based on an analogy with the filtering mass 
defined in the well-studied case of pressure (Gncdin & Hu^ 
ll998l : lGnedinl200d ). In the latter case, the actual gas fraction 
in non-linear, virialized halos, is close to the filtering mass, 
not the Jeans mass, and the filtering mass is affected by the 
value of the Jeans mass at much earlier times. The reason 
the gas fraction is affected by the Jeans mass at early times 
is that the relationship between the gas temperature and 
the gas density is very indirect; the temperature affects the 
pressure gradient, and thus the acceleration, which affects 
the actual position after a delay. In particular, if the tem- 
perature drops suddenly, the gas that was far away from the 
halo center does not instantly fall inside. In the case of cool- 
ing, there are the additional steps from dissociation of H2, 
which changes its abundance, through the process of cool- 
ing which then affects the temperature. The cooling history 
thus affects the distribution of gas. The LW fiux rises expo- 
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nentially fast with time and was very small at early times. 
Thus, the cooling was fast initially, and the gas cooled and 
started to collapse. A sudden late rise in the flux may not 
be able to stop this collapse. 

Using the characteristic value of J21 with a realistically 
large uncertainty should suffice for our main goal of span- 
ning the possible range of the effect of Jlw and Ubc on the 
21-cm background during the X-ray heating era. Specifically, 
we consider two possible feedback strengths which we refer 
to as "weak" and "strong" feedback. Namely, for halos form- 
ing (i.e., virializing) at some time fvir, we adopt the effective 
LW flux J21 in Eq. ([T]) as the LW flux in the same pixel at 
an earlier time tmid, i-e., at the midpoint of the halo forma- 
tion process. In order to obtain a realistically large range 
of uncertainty, with the spherical collapse model in mind 
we either assume that "formation" spans the beginning of 
expansion up to virialization (i.e., t — to t = t^ir, giv- 
ing tmid = |ivir: weak feedback), or just the collapse stage 
starting at turnaround (i.e., t = ^tvir to t — tvir, giving 
tmid = ftvir: stroug f eedback). W e also compare to the lim- 
iting cases (shown in lVisbal et al. ( 2012 )) of no feedback or 
saturated feedback. The latter correponds to assuming that 
star formation is only possible via atomic cooling; in this sce- 
nario, the LW feedback is so efficient that H2 is completely 
dissociated early on and stars form in atomic cooling halos 
(Mcooi > 3 X 10''Mq), as opposed to Mcooi < W'^Mq. For 
reference, we also consider the no feedback case without the 
streaming velocity, in order to assess the importance of the 
velocity effect. For given parameters, at each redshift the 
cosmic mean gas fraction in stars decreases in the different 
cases in the order; no feedback no velocity, no feedback, weak 
feedback, strong feedback, and saturated feedback (where all 
cases except the first include the streaming velocity effect). 

In our hybrid simulations we also incorporate two ele- 
ments of the astrophysics that make our calculations more 
complete and accurate. One aspect is that we include a grad- 
ual low-mass cutoff for star formation, rather than a sharp 
cutoff at Mcooi as we (and others) have previously assumed. 
Since the cooling rate declines smoothly with virial temper- 
ature, a smooth cutoff is expected physically, and indeed 
iMachacek et al.l l|200ll ') found that the fraction of highly- 
cooled, dense gas in their simulated halos is well described 
as being proportional to log(M/Mcooi)- Since this is the gas 
that can participate in star formation, we incorporate this 
by generalizing the star-formation efficiency to include a de- 
pendence on halo mass, /*(M). We assume our standard 
efficiency of /, = 10% for M > Matomic, where Matomic is 
the minimum mass for atomic cooling (~ 3 x IO^Mq but 2- 
dependent). In order for /*(M) to be a continuous function, 
we thus set 

C /. if M ^ Anatomic 

= l f^^-^^LIM^ ffM_i<M<Mato„.c 
I, otherwise. 

(2) 

As shown in Figure [T] (left panel), the standard assump- 
tion of constant /* makes halos with masses near Mcooi 
dominate the cosmic star formation rate, particularly at the 
highest redshifts. Our more realistic model significantly re- 
duces the overall star formation rate (by a factor of 2.0 in the 
example shown a.t z = 19.6) and shifts the peak of the con- 
tribution to star formation to a higher mass (8.7 x Mcooi at 



2 — 19.6). Also shown in the Figure (right panel) is the over- 
all effect of the relative velocity «bc broken down by mass. 
Since the velocity effect on halos is made up of three dis tinct 
effects, with two of them dominant l|Fialkov et al.ll201ll ). the 
dependence on halo mass shows two separate regimes. Near 
the cooling mass (and up to a factor of ^ 2 above it), the 
velocity effect is very strong and also strongly dependent on 
M, mainly due to the boosting of the cooling mass in re- 
gions with a high Ubc. At higher masses, however, the veloc- 
ity effect is weaker and only changes rather slowly with halo 
mass, mainly due to the suppression of the halo abundance. 
A small but non-negligible effect remains even well above 
Matomic. Sluce the velocity effect is strongest at the low-mass 
end (right panel), the shifting of the star formation towards 
higher masses (left panel) reduces somewhat the overall in- 
fluence of the supersonic streaming velocities. Since the LW 
feedback also affects low masses flrst, the log(A/) modula- 
tion delays the LW feedback. 

The second aspect of realistic astrophysics that we in- 
corporate is more directly related to the LW feedback. The 
LW photons emitted by each source are absorbed by hy- 
drogen atoms as soon as they redshift into one of the Ly- 
man lines of the hydrogen atom; along the way, when- 
ever they hit a LW line they may cause a dissociation of 
molecular hydrogen. Some pre vious papers (|Ahn et al.ll2009l : 
iHolzbauer fc Furlanettol l201lh assumed a flat stellar spec- 
trum in the LW region and a flat absorption proflle over 
the LW frequency range. We incorporate the expect ed stel- 
lar sp ectrum of Population III stars from [Barkana fc Loebl 
(|2005h (based on lBromm et al.l l|200ll )'). which varies in the 
LW region typically by a few percent but up to 17%. More 
importantly, we exp licitly include t he fu ll list of 76 rele- 
vant LW lines from iHaiman et al.l l| 19971 ). We summarize 
the results with /lw, the relative effectiveness of causing 
H2 dissociation via stellar radiation. Speciflcally it is the 
ratio between the dissociation rate of molecular hydrogen 
and the naive total stellar flux (i.e., calculated without any 
absorption and integrated over all wavelengths), normalized 
to unity in the limit of zero source-absorber distance. This 
quantity is simply a function of source-absorber distance at 
each redshift under the simplifying assumption of a universe 
at the mean density. This assu mption follows o ur ap proach 
for X-rays (as in 21CMFAST. iMesinger et al.l HoTJ)), and 
should be sufficiently accurate since the strong bias of star- 
forming halos at these high redshifts implies that fluctua- 
tions in star formation (which drive the 21-cm fluctuations) 
are much larger than the fluctuations in the underlying den- 
sity. Thus, give our assumed stellar spectrum we can pre- 
calculate /lw and include this as an effective optical depth 
that is spherically symmetric around each source. Any such 
symmetric effect is easily incorporated within the numeri- 
cal method of 21CMFAST which uses Fourier transforms to 
rapidly perform averages over spherical shells. 

FigureOshows /lw versus the absorber-source distance; 
we parametrize this distance in terms of the absorber-source 
scale-factor ratio R, since /lw versus R is independent of 
redshift. Beyond the max shown R — 1.054 (which corre- 
sponds to 104 comoving Mpc at 2 = 20), /lw immediately 
drops by flve orders of magnitude. The Figure shows that 
LW absorption is poorly approximated as being uniform in 
frequency. In reality, emission from distant sources is ab- 
sorbed more weakly. For example, in one of our main exam- 
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M/Mo M/M, 



Figure 1. Star-formation contribution and effect of velocities versus halo mass (LW feedback not included). Left: The logarithmic 
contribution of each halo mass to the total fraction of gas in stars (i.e., d/55toiiar/'^log(M) averaged over the distribution of I'bc): including 
the log(Af) modulation in Equation JJJl (solid) or with the standard assumption of a fixed efficiency with mass (dashed). We consider 
2 = 13.6 (red), z = 19.6 (green), and z = 25.6 (blue). Right: The ratio of the cosmic mean stellar fraction with v^^^ to the value without 
the velocity effect, i.e., < /stellar (^> f be) >i'bo //stellar (A/, 0). We include Equation l(2j, and consider the same redshifts as in the left 
panel. 
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Figure 2. The relative effectiveness of causing H2 dissociation 
in an absorber at Za due to stellar radiation from a source at 
Zs, shown versus the ratio R = (1 + Zs)/(1 + Za) (solid). For 
comparison we sh ow fmod i ^ commonly used approximation from 
lAhn et al. I 1I2OO9I') (dashed). Both functions are normalized to 
unity at i? = 1. (There was also a 1.4 5% normalization differ - 
ence after we carefully normalized as in lMachacek et al. I l l200j) . 
since we use their results for the LW feedback.) 



pies in the following section (i.e., our strong LW feedback 
case including velocities), assuming a uniform spectrum and 
absorption profile at z = 20 would imply that typically, an 
atom receives 50% of its LW flux from sources out to a dis- 
tance of 18.9 Mpc, 80% from up to 42.2 Mpc, and 90% from 
up to 55.8 Mpc. Our more accurate /lw reduces these num- 
bers to 14.4, 33.0, and 46.2 Mpc, respectively. The accurate 
/lw reduces the overall LW intensity by ~ 20% (thus delay- 
ing the LW feedback), and makes it more short-range (i.e., 
local) and variable. 



3 RESULTS 

3.1 Mean Evolution 

The relative velocity amplifies heating fiuctuations in the 21- 
cm power spectru m, making it possib l e to observe the BAO 
in th e first stars (Visb al et all |2012| : iMcQuinn fc O'Learvl 
l2012h . In this section we consider the effect of the negative 
LW feedback on this exciting observational prospect. Our 
simulation evolves a realistic sample of the un iverse from 
z = 60, roughly when t he first stars turn on l|Naoz et al.l 
I2OO6I : iFialkov et al.|[20Tl! '). We follow a box that is 384 Mpc 
on a side (all distances comoving), with a pixel size of 3 Mpc. 
The X-ray and LW backgrounds in each pixel are made up 
of contributions by stars located within the corresponding 
effective horizons. Since we focus on the era after Lyman- 
a coupling but prior to significant reionization, the 21-cm 
brightness temperature (relative to th e cosmi c micr owave 
background (CMB) temperature TnMR lMadaul |l993) 

rb = 40(l + 5)(l-^)yi±imK, (3) 

where 5 is the gas overdensity and Tk its kinetic tempera- 
ture. 

The negative feedback suppresses star formation on av- 
erage, which leads to a slower rise of the radiative back- 
grounds, and delays various milestones of the star formation 
history. Figure |3] shows the rise of the mean LW flux with 
time. Comparing the two no-feedback cases, we see that the 
streaming velocity has a large overall suppression effect at 
high redshifts, which reaches about an order of magnitude 
at z > 40 but becomes quite small at z < 15. Feedback 
can potentially be very strong at high redshifts (as indi- 
cated by the saturated feedback case), but in practice the 
LW feedback is expected to begin only when the effective 
fiux reaches a level of J21 ~ 10~^; this happens at around 
redshift 30 (weak feedback) or 40 (strong). In both realis- 
tic feedback cases, the LW feedback effectively saturates at 
z - 10. 

We can easily understand why the two feedback cases 
converge with time. Initially, the effective LW flux for star 
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1+z 

Figure 3. The actual LW intensity (solid lines) and the effective 
LW intensity for feedback on star formation (dashed; shown only 
for the two realistic feedback cases). We show the cosmic mean 
intensity (i.e., averaged over our box) versus 1 + z in the follow- 
ing cases: no feedback no vy^^ (purple), and with vy^^: no feedback 
(red), weak feedback (blue), strong feedback (green), and satu- 
rated feedback (black). 

formation [i.e., J21 in Eq. is much higher at a given z for 
the strong feedback case (which assumes a less delayed value, 
closer to the value of J21 at z). The strong resulting feedback 
leads to a slower rise of the actual J21 and thus, eventually, 
also of the effective J21, compared to the weak feedback 
case. Therefore, the effective J21 in the weak feedback case 
gradually catches up with the strong feedback case. Also 
important is that the rate of increase of the flux naturally 
slows with time (i.e., the curves flatten), since star- forming 
halos become less rare (i.e., they correspond to less extreme 
fluctuations in the Gaussian tail of the initial perturbations). 
The weak feedback case effectively looks back to J21 at an 
earlier time, when the rise was faster. 

Figure [3] tracks the rise of the LW flux through several 
milestone values. A reasonable definition of the central red- 
shift zuw of the LW transition is a mean effective intensity of 
J21 =0.1, at which the minimum halo mass for cooling (in 
the absence of streaming velocities) is raised to ~ 2 x 10® 
due to the LW feedback. This is a useful fiducial mass scale, 
roughly intermediate (logarithmically) between the cooling 
masses obtained with no LW flux and with saturated LW 
flux. The central range of the LW feedback transition can 
be defined by the effective LW flux coming within an order 
of magnitude of its central value, so that the minimum Mcooi 
goes from 8 x W^Mq to 5 x IO'^Mq during this period. 

Feedback also slows down the heating of the Universe 
(Figure |4]). For example, the average heating rate at red- 
shift 20 for the weak, strong and saturated feedbacks are 
55.9%, 33.7% and 19.1% of the heating rate with no feed- 
back (all including the streaming velocity). As a result, the 
heating transition is delayed. There are two possible natural 
definitions for this transition, the standard more physical 
definition as the redshift Zh when the mean gas temperature 
equals that of the CMB, and the more observational (or 21- 
cm-centric) definition as the redshift (which we denote zq) at 
which the cosmic mean 21-cm brightness temperature van- 
ishes (Tb) = 0. We consider both definitions, but due to our 
focus on observational predictions, we mostly use zq. 

In our simulation, Zh ~ 17.1 and zq = 16.6 for the no- 




lo 15 20 25 

1+z 



Figure 4. The cosmic heating transition. We show the cosmic 
mean gas kinetic temperature versus 1 -|- 2 in the following cases: 
no feedback no v^,^ (purple), and with vx^^: no feedback (red), 
weak feedback (blue), strong feedback (green), and saturated 
feedback (black). Also shown for comparison is the CMB tem- 
perature (dotted), which crosses the gas temperature at zi^. 

feedback case (no feedback and no velocity gives Zh = 17.7 
and Zq = 17.4), while saturated feedback would delay these 
milestones to Zh = 14.6 and 20 = 14.2. The realistic feedback 
cases are intermediate: Zh = 15.7 and zo = 15.2 for the weak 
feedback case (with a LW transition centered at zlw ~ 19.2, 
and a central range of 2 = 22.0 — 15.2), while Zh = 15.0 
and Zo — 14.6 for strong feedback (with zlw ~ 23.6, and 
a central range of z — 28.3 — 18.1). In every case, the LW 
transition starts very early, and passes through its central 
redshift before the heating transition (with a much bigger 
delay between the two transitions in the strong feedback 
case) . We note that if all the fluctuations were linear, then we 
would find 2h — zo identically. The difference of Az = 0.4 — 
0.5 between them is an example of the effect of non-linear 
fluctuations (plus, in this case, of the non-linear dependence 
of the brightness temperature on the gas temperature). This 
shows that analyses of this era based on linear theory can 
only give rough estimates, and a hybrid simulation like ours 
is necessary in order to properly incorporate the non-linear 
fluctuations in stellar density and other derived quantities. 

3.2 Spatial Fluctuations 

The most interesting 21-cm signature of the first stars is the 
enhanced large-scale fluctuation level due to the supersonic 
streaming velocity. A typical two-dimensional slice of our 
simulated volume is shown in Figure [S] together with the 
spatial distribution of the fluctuations in the density and in 
the relative velocity. A snapshot of the Universe at a fixed 
redshift would look very distinctive in the various feedback 
cases mainly due to the overall delay in the heating due to 
the change in the mean heating rates. It is more instructive 
to compensate for this shift and instead compare each case 
to the others at the same time relative to the individual 
heating redshift (we use zq). 

As seen in Figure [5] the feedback weakens the effect of 
the relative motion, since it boosts the minimal cooling mass 
so that the stars form in heavier halos which are less sensi- 
tive to the relative velocities. Thus, in Figure [5] the 21-cm 
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Figure 5. A two-dimensional slice of our simulated volume. All panels show the same slice, i.e., with the same initial conditions. Top: 
Left: The magnitude of the relative velocity between baryons and dark matter shown in units of the velocity's root-mean-square value. 
Right: The relative fluctuation in density at redshift 20. Bottom: The 21-cm brightness temperature Tj, (relative to the cosmic mean 
in each case) in mK, shown at zq + 3. The cases shown (from left to right) are: no feedback, no v^,^; no feedback; weak feedback; 
strong feedback; and saturated feedback (where the last four include the streaming velocities) . The cosmic mean values (which have been 
subtracted from the maps) are: (T^) = —81 mK, -89 mK, -73 mK, -88 mK and —108 mK, respectively. 



signal with saturated feedback shows the same pattern as in 
the no Vhc case (with no feedback); namely, both of them 
follow the density fluctuations but with a strong enhance- 
ment, where the bias is stronger for the saturated feedback 
(since more massive halos are more strongly biased). The no 
feedback case (with Vhc) shows a strong imprint of the ve- 
locity field along with the influence of density fluctuations; 
e.g., the enhanced Tb at the bottom right is mainly due to 
a density enhancement, while the void at the top (just right 
of center) is mainly due to a large relative velocity (but note 
that the 21-cm maps are the result of a three dimensional 
calculation of the radiation fields, so they cannot be pre- 
cisely matched with two-dimensional slices of the density 
and velocity fields). 

The two realistic feedback cases are intermediate, show- 
ing a clear velocity effect, though not as strong as in the no- 
feedback case. To understand the comparison between the 
weak and strong feedback, we note that the velocities cause 
a very strong suppression of star formation up to a halo mass 
M ~ IO^Mq, but above this critical mass the suppression 
and its M-dependence weaken considerably (Figure [TJ right 
panel). Thus, once the LW feedback passes through its cen- 
tral redshift, the remaining iibc effect changes only slowly 
with M, so that around the time of the heating transition, 
the weak and strong feedback cases show a similar fluctu- 
ation pattern. However, the strong dependence of bias on 
M remains, so that the strong feedback case leads to larger 
fluctuations on all scales. 

These and other features can be seen more clearly and 



quantitatively in the power spectra (Figure [6]). The 21-cm 
fluctuations initially rise with time as the heating becomes 
significant (first in the regions with a high stellar density). 
Eventually, as the heating spreads, the 21-cm fiuctuations 
decline, since the 21-cm intensity becomes independent of 
the gas temperature once the gas is much hotter than the 
CMB (Equation [3]). Thus, the power spectrum reaches its 
maximum height somewhat earlier than zo . The comparison 
among the various feedback cases is complex, since the nega- 
tive LW feedback has several different effects: 1) The lowest- 
mass halos are cut out, reducing the effect of the streaming 
velocity; 2) The higher-mass halos that remain are more 
highly biased; 3) The overall suppression of star formation 
delays the heating and LW transitions to lower redshifts; 
and 4) Since the higher-mass halos that remain correspond 
to rarer fluctuations in the Gaussian tail, their abundance 
changes more rapidly with redshift, making the heating tran- 
sition more rapid (i.e., focused within a narrower redshift 
interval). Thus, at zq + 3 the large-scale {k = 0.05 Mpc~^) 
peak is lower for the realistic feedback cases than it would 
be with no feedback (due to effect #1), and higher for strong 
feedback than for the weak case (due to effect #2). Further 
back in time {zq + 9), weak feedback gives a higher large- 
scale peak than both strong feedback (due to effect #4) and 
no feedback (due to effect #2); at that redshift, saturated 
feedback shows no velocity effect (due to effect #4). 

Lower redshifts offer improved observational prospects, 
due to the lower foreground noise, which makes negative 
feedback advantageous due to effect ^^3, above. We find 
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Figure 6. Power spectra of the 21-cm brightness temperature for no feedback (red), weak feedback (blue), strong feedback (green) and 
saturated feedback (black). Left: z = zq + 3 with (solid) or without (dashed) the relative velocity. Right: Including v^,^, at redshifts 
2 = 20 (dashed), z = 20 + 3 (solid) and 2 = 20 + 9 (dotted). 



that the most promising redshift is z ~ + 3 (Table 1). 
Assuming a first-generation radio telescope arr ay with a 
noise power spectrum th at scales as (1 -I- z)^'^ (|McQuinnl 
I2OO6I : IVisbal et ai]|2012h . the maximal S/N of the large- 
scale {k = 0.05 Mpc~^) peak is 3.24 for weak feedback (at 
z = 18.3) and 3.91 for strong feedback (at z ^ 17.7). For 
compa rison, the no-feedback case considered in IVisbal et al.l 
l|2012l) gave (at z = zo = 2Q) a. S/N of only 2.0 . It would 
be particularly exciting to detect the evolution of the 21-cm 
power spectr um througho ut the heating transition, as we 
suggested in IVisbal et al. I (goij). The S/N remains above 
unity at all z < zq + 7.9 = 23.1 in the case of the weak 
feedback and z < zq + 7.2 = 21.9 for the strong feedback 
(down to 2: = 10 where our simulations end). The streaming 
velocity clearly plays a key role in creating this extended 
observable redshift range, by boosting the large-scale power 
(Figure El). 

We suggested in IVisbal et al.l (|2012h that beyond just 
detecting the power spectrum, it would be particularly re- 
markable to detect the strong BAO signature, since this 
would confirm the major influence of the relative velocity 
and the existence of small (W^Mq) halos. We find that the 
S/N for the large-scale BAO feature of the power spectrum 
is typically ~ 0.5 — 0.7 times that of the large-scale peak 
itself (Table 1). In particular, the BAO S/N also peaks at 
zo + 3, exceeds unity at zq — 0.7 < z < zq + 6.9 (weak) and 
Zq ~ 1.1 < z < Zq + 6.4 (strong) and reaches a maximum 
value of 1.79 (weak) or 2.14 (strong feedback). 

We have assumed here the projected sensitivity of a 
thousand-hour integration ti me with an instrumen t like the 
Murchison Wide-field Array (|Bowman et al.l [20091 ) but de- 
signed to operate in the range of 50-100 MHz. An in- 
strument similarly based on the Low Frequency Array 
iHarker et al.ll2010l ) should improve the S/N by a factor of 
~ 1.5, while a second-generation instrument like the SKA 
or a 5000-anten na MWA should improve it by at least a 
factor of 3 or 4 (|McQuinnll200^ : IVislpal et al.ll2012l ). Thus, 
future instruments may be able to probe even earlier times, 
including the central stages of the LW feedback. 



4 CONCLUSIONS 

We have presented new predictions for the signature of the 
first stars in the heating fluctuations of the 21-cm bright- 
ness temperature. We ran hybrid simulations that allow 
us to predict the large-scale observable 21-cm signature 
while accounting on small-scales for various effects on star- 
formation investigated by previous analytical models and 
numerical simulations. In particular, we incorporated for the 
first time the Lyman- Werner feedback on star formation, 
calculated self-consistently including the effect of the super- 
sonic streaming velocity. A three-dimensional calculation of 
the Lyman- Werner and X-ray backgrounds allowed us to 
calculate the heating history of the gas and the resulting 
21-cm intensity maps. 

We have focused on the negative LW feedback, which 
begins at z ~ 30 — 40 but strengthens very gradually, passing 
its central point at zlw ~ 19 — 24 and saturating only at z ~ 
10. The heating transition is centered at zo ^ 15 (including 
a delay of Az ~ 1.5 — 2 due to the feedback), when the LW 
transition is well advanced but still far from saturated. The 
large-scale 21-cm power spectrum is potentially observable 
over a broad redshift range of z ~ 10 — 22 or 23. The best 
prospects are at z ^ 18, when the large-scale peak reaches 
a signal-to-noise ratio (for a projected first-generation radio 
telescope array) of 3.2 (for our weak feedback case) or 3.9 
(for strong feedback). At this redshift, the BAO signature 
(which marks the velocity effect and the presence of 10^ M© 
halos) should also be observable with a S/N~ 2. The BAOs 
should be observable over a broad redshift range of Az ~ 
7.5 . 

These numbers are obtained with our standard set of ex- 
pected astrophysical parameters, but they may shift around 
a bit depending on the precise properties of the early stars 
and their remnants. We hope these findings will stimulate 
additional numerical simulations of the complex radiative 
feedback at z ~ 10 — 30, as well as future observational 
efforts in 21-cm cosmology directed at the epoch prior to 
reionization. 
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Z - ZQ 


no v^^ 


5Tb (A: = 0.05 Mpc" 
no fbk weak 


■1\ q/l\r 
strong 




no tijjc 


BAO, S/N 
no fbk weak 


strong 


Sclt 


-3 


1.07 


1.24 


1.60 


1.73 


1.84 


0.45 


0.58 


0.72 


0.76 


0.79 





1.68 


2.33 


2.35 


2.69 


3.09 


0.70 


1.26 


1.16 


1.30 


1.31 


3 


2.26 


3.59 


3.24 


3.91 


4.74 


0.91 


2.18 


1.79 


2.14 


2.00 


6 


1.02 


1.75 


2.08 


1.89 


1.34 


0.37 


1.17 


1.30 


1.18 


0.54 


9 


0.086 


0.33 


0.56 


0.34 


0.31 


0.051 


0.23 


0.41 


0.25 


0.14 


12 


0.18 


0.23 


0.24 


0.27 


0.34 


0.083 


0.099 


0.11 


0.12 


0.15 



Table 1. The signal to noise ratio S/N (i.e., the square root of the ratio between the power spectra of the signal and noise), for a 
projected first-generation radio array. We show the S/N of the large-scale peak at the wavenumber k = 0.05 Mpc~^ (Left), and of the 
BAO component (Right), at various redshifts, for five cases: no feedback no ti^c- no feedback (with ?Jbc)i weak, strong and saturated 
feedback. The BAO S/N is defined as the square root of the difference between the peak at fc = 0.05 Mpc~^ and the trough at A; = 0.07 
Mpc~^, each measured with respect to the non-BAO power spectrum (i.e., the power spectrum smoothed out using a quartic fit), and 
each normalized by the noise power spectrum at the same k at the corresponding redshift. 
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